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INTRODUCTION 

This  report  is  divided  into  three  sections.  In  the  first  chapter, 
the  accuracy  of  coupled  mode  theory  as  applied  to  volume,  dielectric, 
reflection  gratings  which  have  unslanted  fringes  and  either  periodic  or 
aperiodic  refractive  index  profiles  is  investigated.  The  accuracy  is 
evaluated  by  comparing  coupled  mode  theory  results  with  those  obtained 
by  Abele's  exact  multilayer  theory.  For  the  examples  considered,  coupled 
mode  theory  gave  accurate  results  only  for  gratings  having  periodic 
refractive  index  profiles. 

In  Section  2,  a  method  for  designing  volume,  dielectric,  reflection 
gratings  having  unslanted  fringes  is  developed.  The  technique  is  ap¬ 
plicable  to  two  types  of  problems.  The  first  type  is  one  in  which  the  re¬ 
flectance,  R,  vs.  angle  of  incidence  and  wavelength  is  specified,  and  the 
second  type  is  one  in  which  the  amplitude  reflection  and  transmission  co¬ 
efficients,  r  and  t,  respectivly,  vs.  angle  of  incidence  and  wavelength  are 
specified.  In  both  cases,  the  technique  determines  (approximately)  the  cor¬ 
responding  one  dimenstional  refractive  index  profile,  n(z).  The  synthesis 
method  is  illustrated  by  two  examples,  and  for  these  examples,  the  method 
is  seen  to  be  reasonably  accurate.  In  the  appendix  to  Section  2,  the  re¬ 
flection  grating  design  approach  is  extended  to  multilayer  dielectric  filters. 

In  Section  3,  a  formal  mathematical  analogy  between  reflection  grat¬ 
ings  and  corrugated  waveguide  filters  (CWF)  is  demonstrated.  The  possi¬ 
bility  of  designing  CWF,  using  the  iterative  Fourier  transform  technique 
developed  in  Section  2,  is  explored.  It  is  emphasized  that  the  mathe¬ 
matics  developed  for  CWF  is  not  rigorous.  Consequently,  predicted  re¬ 
sults  should  be  experimentally  verified  and  attempts  should  be  made  to 
develop  a  rigorous  mathematical  approach. 
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Section  1 


ACCURACY  OF  COUPLED  MODE  THEORY 


In  this  section,  we  examine  the  accuracy  of 
coupled  mode  theory  as  applied  to  volume,  dielectric, 
reflection  gratings  which  have  unslanted  fringes  and 
either  periodic  or  aperiodic  refractive  index  pro¬ 
files.  The  accuracy  is  determined  by  comparing  coupled 
mode  theory  results  with  those  obtained  from  Abelfe's 
exact  multilayer  theory.  For  the  examples  considered, 
coupled  mode  theory  gave  accurate  results  only  for 
gratings  having  periodic  refractive  index  profiles. 
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1-1.  INTRODUCTION 

The  study  of  wave  propagation  through  periodic  media  dates  back 
almost  one  hundred  years,  and  has  had  a  significant  impact  in  a  number 
of  areas  including  the  theory  of  electron  energy  bands  in  crystals, 
microwave  devices,  diffraction  gratings,  and  thin  film  optical  devices. ^ 
Analytic  solutions  to  these  wave  propagation  problems  were  generally 
based  on  Floquet's  Theorem. ^  During  the  1960's,  considerable  effort 
was  directed  toward  finding  exact  solutions  of  optical  propagation 
through  sinusoidally  periodic  volume  gratings.  In  particular,  C.B. 
Burckhardt  solved  the  problem  of  diffraction  of  a  plane  wave  by  a  planar, 
sinusoidal,  dielectric  transmission  grating  having  unslanted  fringes. ^3 
Unfortunately,  the  solutions  obtained  were  not  in  closed  form  and  in¬ 
volved  approximate  nunerical  techniques  for  their  evaluation.  These 
numerical  techniques  were  often  nontrivial  to  implement.  The  numerical 
difficulties  encountered  together  with  the  lack  of  insight  gained  from 
these  numerical  techniques  was  an  impetus  for  developing  a  simpler  approach. 

f4  5i 

A  natural  candidate  was  coupled  mode  theory. L  ’  J  This  theory  had 
already  been  used  quite  sucessfully  to  study  microwave  devices. ^ 

In  many  instances,  one  is  interested  in  a  structure  into  which  a  small 
perturbation  has  been  introduced.  The  modes  which  propagate  in  the  non- 
perturbed  structure  are  often  easily  computed.  The  perturbation  causes 
some  of  these  modes  to  interact  with  one  another.  Generally  the  structure 
is  excited  in  such  a  manner  that  a  single  mode  is  launched  into  the 
structure,  and  the  perturbation  causes  this  mode  to  excite  and  exchange 
energy  with  other  modes  as  it  propagates.  In  coupled  mode  theory 
approximate  equations  are  derived  to  describe  this  coupling  between 
modes . 

Kogelnik  used  coupled  mode  theory  to  derive  closed  form  expressions 
for  plane  wave  Bragg  diffraction  from  thick,  planar,  sinusoidal  hologram 
gratings.  ^3  The  gratings  could  have  arbitrary  fringe  slant  and  operate 
in  either  the  reflection  or  transmission  mode.  Kogelnik's  results  were 
later  extended  to  periodic  gratings  having  arbitrary  shape  (i.e.  not 
necessarily  sinusoidal ). ^ 
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Unlike  Floquet's  Theorem,  the  coupled  mode  approach  does  not  re¬ 
quire  that  the  perturbation  be  periodic.  There  has  been  some  interest 
in  producing  corrugated  waveguide  devices  for  filtering  which  have 
nonperiodic  corrugations.  The  characteristics  of  these  structures 
have  been  analyzed  using  coupled  mode  theory. appears  quite 
difficult  to  determine  bounds  on  the  accuracy  of  the  coupled  mode 
approach.  We  note,  however,  that  Bragg  diffraction  from  periodic  or  non¬ 
periodic  reflection  holograms  having  zero  fringe  slant  can  be  solved 
exactly  using  a  simple  multilayer  matrix  technique  developed  by  Abelfc 
in  the  1950's.  J  Thus,  these  reflection  grating  structures  can  be 

used  as  a  check  on  the  accuracy  of  the  coupled  mode  approach.  A  re¬ 
flection  hologram  of  the  type  discussed  above  is  analogous  to  corrugated 
waveguide  devices  and  consequently  takes  on  added  importance. ^ 

Sakuda, Stoll  and  Yariv  have  compared  coupl ed  wave  theory  results 
and  exact  results  (derived  via  Floquet's  Theorem)  for  sinusoidally  cor¬ 
rugated  waveguides.^-1 6-1 Their  analysis  indicated  that  the  two 
approaches  agree  reasonably  well.  Moharam  and  Gaylord  also 
ran  comparisons  for  unslanted,  sinusoidal,  reflection  and  transmission 
holograms.  For  reflection  holograms  having  unslanted  fringes,  their 
analysis  indicated  that  the  two  approaches  give  results  which  are  in 
excellent  agreement.  For  slanted  fringes  they  show  substantial  dis¬ 
agreement,  but  they  have  assumed  very  high  refractive  index  modulations 
(An/n  as  high  as  0.29). 

In  this  section  we  present  four  examples  in  which  we  compare  results  of 
coupled  wave  theory  with  those  of  Abeles'  exact  multilayer  theory  for  re¬ 
flection  gratings  having  unslanted  fringes.  The  examples  show  that  there 
is  substantial  disagreement  between  the  two  approaches  for  gratings  hav¬ 
ing  nonperiodic  refractive  index  profiles. 

1-2.  PROBLEM  STATEMENT 

Consider  a  pure  dielectric  material  which  has  a  refractive  index 
variation,  n(z),  only  along  the  z-axis  as  shown  in  Figure  1-1. 
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FIGURE  1-1.  DIELECTRIC  MATERIAL  HAVING  A  REFRACTIVE 
INDEX  VARIATION  IN  THE  Z  DIRECTION.  Also 
shows  are  incident,  reflected  and  transmitted  beams. 
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The  refractive  index  variation  is  written  as 


n(z)  =  nQ  +  An(z)sin(27rf(z)z) 


(1-1) 


where 


An(z)«  nQ 


(1-2) 


The  material  is  assumed  to  be  of  infinite  extent  in  both  the  x  and  y 
dimensions.  Regions  c  and  S  shown  in  Figure  1-1  have  refractive  indices 
nc  =  n(o)  and  n$  =  n(D)  respectively.  An  s  polarized  (electric  field 
along  the  y  axis)  plane  wave  of  wavelength  \  is  incident  upon  the  material. 
Some  of  this  incident  field  is  reflected,  and  the  remainder  is  trans¬ 
mitted.  The  incident  electric  field  amplitude  is  designated  ER,  while 
the  reflected  and  transmitted  electric  field  amplitudes  are  designated 
EL  and  Er  resPect1'vely-  We  wish  to  determine  the  reflection  and  trans¬ 
mission  coefficients  r  and  t  defined  by 


( 1  -3a ) 
( 1  -3b ) 


The  general  procedure  for  determining  r  and  t  is  to  solve  Maxwell's 
Equations  in  the  three  regions  shown  in  Figure  1-1,  and  then  match  the 
tangential  electric  field  and  magnetic  field  components  of  these  solu¬ 
tions  at  the  boundaries  separating  regions  C  and  G  and  regions  G  and  S. 

We  proceed  as  follows:  The  wave  equation  in  all  three  regions 
can  be  written  as 

2 

V2Ey(x,z)  +  (£j)  n2(z)Ey(x,z)  =  0  (1-4) 

where  £y  is  the  electric  field  component  along  the  y  axis.  Since 
An(z)  «  nQ,  we  can  write 

n2(z)  a  n2  +  2nQAn(z)  sin  2nf(z)z  (1-5) 
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In  region  G,  Eq.  (1-4)  has  a  solution  of  the  form 
Ey(x,z)  =  A(z)exp(-i(px  x  +  p2  z))  +  B(z)exp(-i (px  x  -  pz  z ) ) 
where 

^  2tt 

px  =  -nQ  sine 

A  2tt  - 

o2  '  T  "o  cose 

Substituting  Eqs.  (1-5)  and  (1-6)  into  Eq.  (1-4)  yields 

A"(z)exp(-1(pxx  +  pzz))  -  2ipz  A' (z)exp(-i (Pj{x  +  pzz)) 

-  pz  A(z)exp(-i(pxx  +  pzz))  -  px  A(z)exp(-i(pxx  +  Pzz)) 

+  B"(z)exp(-i(pxx-  pzz))  +  2ipz  B' (z)exp(-i (pxx  -  Pzz)) 

-  Pz  B(z)exp(-i(Pxx  -  p2z))  -  px  B(z)exp(-i(pxx  -  Pzz)) 

2  2 
+  (~  nQ)  A(z)exp(-i(pxx  +  Pzz))  +  (H  nQ)  B(z)exp(-i (pxx  - 

-  i (^fnp An(z)  A(z)exp(-i(pxx  +  pzz)  +  i2irf(z)z) 

2 

+  i(^"0)  &n(z)  A(z)exp(-i(pxx  +pzz)  -  i2irf(z)z) 

-  i(^fnQ^n(z)  8(z)exp(-i  (p  xx  -pzz)  +  i2irf(z)z) 

+  i(%2n0An(z)  B(z)exp(-i(pxx  -pzz)  -  i2irf(z)z)  =  0 
Using  Eqs.  (l-7a)  and  (l-7b),  Eq.  (1-8)  reduces  to 

A"(z)exp(-i(pxx  +  p2z)  -  21 pz  A' (z)exp(-i (pxx  +  pzz)) 

+  B"(z)exp(-i(Pxx  -  p2z))  +  21pz  B ' (z)exp(-i (pxx  -  pzz)) 


(1-6) 

( 1  -  7a ) 
(l-7b) 

P2*)> 

(1-8) 
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"  1 (^f)2n0^n(2)  A(z)exp(-i(pxx  -  pzz)  +  iA3(z)z) 

+  i(^ncAn(^  A(z)exp(-i(pxx  +  3pzz)  -  IaB(z)z) 

'  1  (^)2n0An(z)  B(z)exp(-i(pxx  -  3pzz)  +  iAfc>(z)z) 

+  i(^-)  noAn(z)  B(z)exp(-i(pxx  +  pzz)  -  iA3(z)z)  ■  0  (1-9) 

where 

AB(z)  -  2nf  (z)  -  2pz  (1-10) 

It  is  important  to  note  that  no  approximations  have  been  made  in 
the  derivation  leading  up  to  Eq.  (1-10). 

1-3.  COUPLED  MODE  THEORY 

It  is  common  in  coupled  mode  theory  to  derive  from  Eq.  (1-9)  two 
coupled  mode  differential  equations  when  Ap(z)  is  small.  First,  it 
is  assumed  that  only  termo  in  Eq.  (1-9)  which  have  approximately  the 
same  phase  factors  will  significantly  interact  with  one  another. 

Second,  it  is  assumed  that  A(z)  and  B(z)  change  phase  slowly  as  com¬ 
pared  to  exp(-ipzz).  With  these  assumptions  Eq.  (1-9)  becomes 

A" (z)  -  2ipz  A*  (z)+ i  {^-)  nQAn(z)exp(-iA3(z)z)B(z)  =  0  (1-lla) 

p  2 

B"(z)  +  2ip z  B'(z)-i(^)  n0An(z)exp(iA3(z))A(z)  =  0  (1-llb) 

Eds-  (1-lla)  and  (1-llb)  are  further  simplified  by  making  the  assumption 
that 


| A" (z ) | << |p z  A'(z)|  and  jB"(z)|  «  |pz  B’(z)|  (1-12) 


With  the  above  approximation,  Eqs.  (1-lla)  and  (1-llb)  become 
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1  o  2 

A'(z)  =  (=^)  n0An(z)exp(-iAB(z)z)B(z)  ( 1  - 1 3a ) 

1  2 

B'(z)  =  2p~  nQAn(z)exp(iA8(z)z)A(z)  ( 1  -1 3b) 

Eqs.  (l-13a)  and  (1-1 3b)  are  the  first  order  coupled  mode  equations 
which  commonly  appear  in  the  literature.  Since  the  refractive 
index  is  continuous  across  the  boundaries  between  regions  C  and  G  and  re 
gions  G  and  S,  Maxwell's  boundary  conditions  require 


> 

o 

II 

m 

23 

( 1  - 1 4a ) 

6(0)  •  El 

(1-1 4b ) 

A(D)  -  \ 

(1-1 4c ) 

6(D)  •  0 

(l-14d) 

(1-1 4e ) 

(1-1 4f ) 

1-4.  ABELES'  EXACT  MULTILAYER  THEORY 

In  the  1 950 ' s  Abeles  developed  a  convenient  numerical  matrix 

technique  for  solving  Eq.(l-9)P^  First,  the  structure  shown  in 

Figure  1-1  is  divided  into  thin  slabs  each  parallel  to  the  x-y  plane. 

The  slabs  are  chosen  to  be  sufficiently  thin  so  that  the  refractive 

index  n(z)  can  be  assumed  constant  within  each  thin  slab.  Thus,  a 

solution  of  the  wave  equation  within  the  mth  thin  slab,  which  is  assumed 

to  have  a  constant  refractive  Index  of  n  ,  is 

m 

VX,2)  ■  +  l'-'5) 
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where 


n  -  n 

pxm  x  nm 


4  2tt  _ 
pzm  X  nm 


sin  em 
cos  em 


( 1  -1 6a ) 
(1 -1 6b) 


For  m  =  0  and  m  =  N+l  Eq.  (1-15)  is  also  a  solution  of  the  wave  equation 
in  the  half  space  regions  C  and  S  respectively. 

By  requiring  the  tangential  components  of  the  electric  and  magnetic 
fields  to  be  continuous  across  the  boundaries  between  the  slabs, 
it  is  straightforward  to  show  that 


Bm 


where 


V“]  v  U 
m  vm+l  m+1 


/\n+l 

Bm+1 


Vlsin  0m+l  *  Vinem 


n  cos  g 

Ym  =  - z -  is  the  characteristic  impedance 


of  free  space) 


<i)m  =  n„  dm  cos  e„  (d„  is  the  thickness  of  the 
Tm  xmm  m  m  .  ... 

nth  thin  slab) 


v_  * 

m 


u  * 

m 


1  1 


Y  -Y 
L  m  mJ 


exp(i0m) 


mv 


exp<-i*m)J 


(1-17) 

0-18) 

(1-19) 

(1-20) 

(1-21) 

(1-22) 


Now  with  the  following  definitions 
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-  vuv 

m  m  m 


-1 


it  immediately  follows  that 


S 


AN+1 

BN+1 


] 


(1-23) 

(1-24) 


0-25) 


(1-26) 


Furthermore  in  Fig.  1-1, 


Therefore, 


=  0 


(1-27) 


(l-28a) 


( 1 - 28b ) 


1-5.  EXAMPLES 

We  have  calculated  r  vs  X  for  five  different  refractive  index 
profiles  using  both  the  exact  multilayer  technique  and  the  coupled 
mode  approach.  The  coupled  mode  equations  (1-I3a)  and  (l-13b)  were  I 

solved  numerically  using  Runge-Kutta.  The  results  are  shown  in 

<> 

: 
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Figure  1-2.  The  multilayer  technique  gives  exact  results  of  course,  only 
if  the  refractive  index  is  constant  within  each  of  the  slabs.  This  will 
be  the  case,  in  the  limit,  as  each  of  the  slab  thicknesses  go  to  zero. 

The  solid  curves  given  in  Figure  1-2  were  obtained  by  decreasing  all  the 
slab  thicknesses  until  further  reductions  had  no  appreciable  affects  upon 
the  multilayer  theory  predictions.  Excellent  agreement  is  seen  to  exist 
for  the  periodic  grating  but  the  agreement  is  poor  for  the  nonperiodic 
gratings. 

Several  assumptions  have  been  made  in  the  derivation  of  the  first 

order  coupled  mode  equations  (l-13a)  and  (1-1 3b ) .  Gaylord  helieves  that 
the  assumption  given  by  Eq.  (1-12)  may  contribute  the  most  error.  [18] 
Kogelnik  has  shown  that  this  assumption  is  justified  for  sinusoidal 
gratings,  but  we  feel  that  his  argument  is  inlorrect.  Our  reason¬ 
ing  goes  as  follows: 

For  a  sinusoidal  grating  at  the  Bragg  angle  aB(z)  =  0  and  An(z)  is 
a  constant.  Thus,  Eqs.  (l-13a)  and  (1 -1 3b)  may  be  rewritten  as 


A'(z)  =  KB(z)  (1-2  9a ) 

B'(z)  =  KA(z)  ( 1 -29b ) 


where  K  is  a  constant.  Combining  Eqs.  ( 1  -29a )  and  (1-2%)  yields 

A"  (z)  =  K2A(z)  (1-30) 

Now  for  a  reflection  grating  with  unslanted  fringes  B(D)  =  0  and 
A  (D)  f  0.  Thus  Eq.  (l-29a)  implies 

A'(D)  =  0  (1-31) 

while  Eq.  (1-30)  implies 

A" ( D)  t  0  0-32) 

Therefore, 

j  A"  ( D)  |  >|p2A'(D)|  0-33) 

which  contradicts  the  assumption  of  Eq.  (1-12)  forz«=D. 
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1-6.  CONCLUSIONS 

We  have  shown  by  example  that  coupled  mode  theory  can  give  in¬ 
accurate  results  under  certain  conditions.  Furthermore,  one  of  the 
common  assumptions  used  with  coupled  mode  theory,  i.e.,  neglecting 
second  derivatives,  has  been  shown  to  lack  justification.  It  ap¬ 
pears  that  great  care  must  be  taken  when  attempting  to  use  coupled 
mode  theory  for  analyzing  optical  propagation  through  gratings  and 
corrugated  waveguide  filters. 
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Section  2 

AN  ITERATIVE  FOURIER  TRANSFORM  DESIGN  TECHNIQUE 

In  this  section  an  iterative  technique  for  designing  reflection 
gratings  having  arbitrary  reflectance  vs.  wavelength  (or  angle  of 
incidence)  characteristics  is  presented.  The  technique  is  based 
on  Fourier  transforms  and  is  similar  to  the  multilayer  dielectric 
stack  design  approaches  developed  earlier  by  Delano  and  Sossi. 
Reflection  grating  design,  using  the  iterative  technique,  is  illus¬ 
trated  with  two  examples.  In  addition,  the  reflection  grating  design 
technique  is  shown  in  Appendix  A,  to  be  applicable,  with  slight 
modification,  to  the  multilayer  dielectric  stack  design  problem. 


2.1  INTRODUCTION 

Two  approaches  for  determining  the  reflectivity  characteristics 
of  periodic  and  aperiodic  reflection  gratings,  first  order  coupled 
mode  theory  and  Abeles'  multilayer  theory,  were  presented  in  section 
one.  The  two  approaches  were  seen  to  give  identical  results  for 
sinusoidal  gratings  but  substantially  different  results  for  some 
aperiodic  gratings.  Recall  that  Abeles'  multilayer  theory  is 
known  to  be  exact,  whereas  first  order  coupled  mode  theory  is 
approximate.  A  formal  mathematical  analogy  between  corrugated 
waveguide  filters  (CWF)  and  reflection  gratings  will  be  demonstrated 
in  chapter  three,  thus  raising  the  possibility  of  applying  reflection 
grating  design  and  analysis  techniques  to  corrugated  waveguide  filters. 

It  is  the  aim  of  this  chapter  to  describe  an  analytic  approach  for 
designing  reflection  grating  filters.  The  approach  allows  one  to  deter¬ 
mine  the  refractive  index  profile,  n(z),  which  results  in  some  user 
specified  grating  reflectivity  vs,  wavelength  characteristic.  The 
reflection  gratings  to  be  considered  will  be  structurally  identical 
to  those  described  in  section  one. 

Analytic  design  techniques  for  reflection  gratings  and  cor¬ 
rugated  waveguide  filters  have  been  under  investigation  during  the 
last  ten  years.  Kogelnik  has  computed  the  characteristics  of  some 
aperiodic  CWF,  but  has  not  indicated  any  approach  for  designing 
filters  which  would  have  specified  reflectivity  characteristics.^-1 

Snubert  has  demonstrated  that  step  variations  in  the  periodic 
corrugations  of  a  distributed  feedback  laser  can  result  in  lower 
threshold  gains,  greater  higher  order  mode  rejection,  and  unidirec¬ 
tional  operation.  Matsuhara,  et.  al.  have  even  proposed  a 
general  technique  for  CWF  synthesis. ^3  Kogelnik's  analysis 
approach  and  Matsuhara's  analysis  and  design  technique  are  suspect, 
since  they  depend  on  first  order  coupled  mode  theory.  Furthermore, 
even  within  the  framework  of  first  order  coupled  mode  theory, 

Matsuhara's  approach  requires  one  to  make  several  questionable  assumptions. 
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Design  techniques  for  multilayer  dielectric  filters  have  been, 
and  are  still  actively  being,  investigated,  The  reflection 

grating  design  approach,  to  be  described  below,  relies  heavily 
on  multilayer  dielectric  filter  theory.  It  is  quite  similar  to  some 
of  the  multilayer  dielectric  filter  Fourier  Transform  techniques 
developed  by  DelanoP'^Sossf  ,  Pegis,^^  and  Knittl^.  Refer¬ 
ence  8  provides  an  excellent  overview  of  multilayer  dielectric 
filter  Fourier  synthesis  techniques.  Greenewalt,  et.  al .  has 
also  developed  a  Fourier  synthesis  technique.  In  Appendix  A, our 
reflection  grating  design  approach  is  extended  to  multilayer  di¬ 
electric  filters. 

2.2  REFLECTION  COEFFICIENT  DIFFERENTIAL  EQUATION 

We  start  by  deriving  a  differential  equation  for  the  reflection 
coefficient.  Consider  the  pure  dielectric  shown  in  Figure  1-1  which 
has  a  refractive  index  variation  n(z).  The  dielectric  is  assumed 
to  be  of  infinite  extent  in  both  the  x  and  y  dimensions.  Regions 
c  and  s  have  refractive  indices  n£  =  n (0)  and  n$=  n(D),  respectively. 
An  s  polarized  (electric  field  along  the  y  axis)  plane  wave  of 
wavelength  X  is  incident  at  angle  of  6Q  upon  the  dielectric.  Some 
of  the  incident  field  is  reflected,  and  the  remainder  is  transmitted. 

The  incident  electric  field  amplitude  is  designated  E^,  while  the 
reflected  and  transmitted  electric  field  amplitudes  are  designated 
E^  and  E^,  respectively.  The  reflection  and  transmission  coefficients 
are  defined  by  Equations  ( 1 -3a )  and(l-3b). 

In  all  three  regions,  the  total  electric  (E)  and  total  magnetic  (H) 
field  components  can  be  written  as 

Ey(x,z)  =  U(z)exp(-ipxx) 

Hx(x,z)  =  V(z)exp(-ipxx) 

Hz(x,z)  =  W(z)exp(-ipxx) 
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where 

Px  =  —  n(0)sin  eQ 

(2-ld) 

From  two  of  Maxwell's  Equations 

and 

yxF  =  -  iupfzjff 

(2-2a) 

where 

Vxff  =  T  +  iut(z)F 

(2-2b) 

Ott 

w  =  =-  c  (c  is  the  speed  of  light  in  vacuum) 

A 

y(z)  is  the  permeability 

e(z)  is  the  permittivity 

it  follows 

that 

3E 

~  i  <4J(z)Hx 

(2-3a) 

3E 

^  =  -  iwu(z)Hz 

(2-3b) 

“9T  "  =  iu£(z)Ey 

( 2 -3c ) 

Substituting  Equations  (2-la)  -  (2-lc)  into  Equations  (2-3a) 

-  (2-3c) 

yields 

U'(z)  =  iu)y(z)V(z) 

(2-4a) 

_1PxU(z)  =  -  1u*i(z)W(z) 

(2-4b) 

V* (z)  +  i Pxw(z)  =  iwe(z)U(z) 

(2-4c) 

Equations 

(2-4b)  and  (2-4c)  can  be  combined  to  give 

2 

v'(z)  =  1("  lllz7+ 

(2-5) 
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The  permeability  and  permittivity  of  free  space  will  be  designated 
P0  and  £0>  respectively.  The  quantities  e(z),  Y(z)  and  0(z)  are 
defined  below 


COS  e(z)  =  1  -  sin2  0O 

(2-6a) 

x  A  V2 

y(z)  =  (p(z)e(z)/p0e0)  cos  6(z) 

( 2 -6b) 

0(z)  =j\i{z)/e(z)/< cos  9(z) 

(2-6c) 

„  ^  2ir 

K  _  T 

(2-6d) 

Using  the  above  four  definitions.  Equations  (2-4a)  and  (2-5)  can  be 
rewritten  as 

U'(z)  =  i Ky(z)g(z)  V(z)  ( 2 -7a ) 

V*  (z)  =  iK(y(z)/e(z)).U(z)  (2-7b) 

The  input  impedance,  Q(z)  is  defined  by 

QU>  -  yffy  (2-8) 

If  follows  from  Equations  (2-7a)  and  (2-7b)  that 

Q'(z)  =  iK  y(z)0(z)  -  iK(y(z)/g(z))Q2(z)  (2-9) 

Assume  that  the  section  of  the  dielectric  to  the  left  of  the 
plane  z  =  zQ  (0  <  zQ  _<  D)  is  replaced  by  a  homogeneous  dielectric 
of  refractive  index  n(zQ).  Further  suppose  that  an  s  polarized 
plane  wave,  having  x  propagation  constant,  px,  and  electric  field  strength 
Er,  Is  incident  from  the  left  onto  the  z  =  zQ  plane.  A  portion  of  the  in¬ 
cident  plane  wave  is  reflected  and  the  remainder  is  transmitted.  From  Eq. 
(2- 3a)  and  (2-3b)  the  incident  (1)  and  reflected  (r)  plane  wave  fields  can 
be  written  as 
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Eiy(x,z)  =  ERexp(-i(Pxx+  p2(zQ)z)  (2-10a) 

Hix(x-2)  ■  -  'SiTIJ  er  ex|)(-1(V  +  c2I2o>z» 

0  (2-106) 

Hl2(x>2)  =  -Sinn-  er  exp(-i(»xx  +  p2(zo>2) 

0  (2-10c) 


Ery(x,z)  =  ELexp(-i(pxx  -  P2(zQ)z))  (2-10d) 

Pz(zJ 

Hrx(x,z)  =  ~hTCz~')EL  exP("i(pxx  ‘  Pz(zo)z))(2“10e) 

P„ 

Hrz(x>2)  =  m~TD  EL  ®V(-1(Pxx  -  pz(z0>z)> 

0  (2-1 Of ) 

where  _ _ 

pz(z)  -  +V(^)n2(z0)  -  P2  (2-10g) 


Continuity  of  the  tangential  electric  and  magnetic  field  components 
across  the  z  =  zQ  plane  requires  that 


Eiy(x,zQ)  +  Ery(x,zo)  =  U(z0)exp(-ipxx)  (2-lla) 


and 

Hix(x,z0)  +Hrx^x,zo^  =  v(zo}  exPHPxx)  (2-1 lb) 


Combining  Equations  (2-10a)  through  (2-1  lb)  yields 


1  +  El/Er  i*i(z0) 
+  PZ(Zq) 


U(zq) 

V(zo) 


-  Q(zo> 


(2-12) 
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But  by  Equations  (2-ld),  (2-10g),  (2-6a)  and  (2-6c) 


wm(z0) 


(2-13) 


Using  Equations  (2-13)  and  the  definition  of  the  reflection  coefficient, 
r,  given  by  Equation  (1 -3a) ,  Equation  (2-12)  can  be  rewritten  as 


1  +  r(zJ 

:  i  +  ho  3Uo)  =  q(Zq) 


(2-14) 


Finally  substituting  Equation  (2-14)  into  Equation  (2-9)  results 
in  the  following  differential  equation  for  r(zQ) 

r'(z0)  =  i2Ky(z0)r(z0)  -  "  r2(z0))  (2'16) 

Equation  (2-15)  was  originally  given  by  Walker  and  Wax^4-'  and  Kofink^-15^, 
When  | r(z  ) |2  «  1,  then  by  Equation  (2-15) 


B'(zo} 


r'(zQ)  -  i2KY(z0)r(z0)  - 


(2-16) 


The  solution,  r(z),  of  Equation  (2-16)  with  equality  assumed  is 

z  z(  P 

r(z)exp(-  i  2K  yY(^)d^)  -  r(o)  =  -  j  exp[-  2iK  J  YU)dc]dp 

0  °(2-17) 

Imposing  the  boundary  condition  r(D)  =  0,  it  immediately  follows  that 

0  P 

r{0)  =  f  exPt-  i2K  /(?)dc3dP  (2-18) 
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The  refractive  index  n(z)  is  given  by 

n(z)  (2-19) 

"  M0  0 

We  restrict  our  analysis  to  nonmagnetic  dielectrics,  and  therefore. 


y(z)  «  y0  (2-20) 

The  definitions  given  by  Equations  (2-6a)  -  (2-6c)  can  now  be 
rewritten  as 


From  Equations  (2- 21a) 


cos2e(z)=  1  -  sin2e 

n2(z)  0 

( 2-21  a ) 

Y(z)  =  n(z)cos  6(z) 

(2-21 b ) 

^  -  feasm 

(2-2 1 c ) 

and  (2-21c) 

1  dg(z)  _  n'(z) 


n'(z) 

T, 


n(z)cos2  0(z)  n(z)[  t  .  Ao )_  s.n2  77 


nr(z) 


f  y(c)dc  =/n( c)n  -  sin2  eo]  '  dc 

*^o  Jo  n  {$) 


(2-22a) 

(2-22b) 


o  •'o 

But  by  Equations  (1-1)  and(l-2) 


(2-23) 
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and  so 


1 

eTzT 


n‘(z) 
nocos  90 


z 

(?)cos  0O d?  *  no(coseo)z 


(2 -24a) 


(2 -24b) 


Finally,  combining  Equations  (2-18)  with  (2-24a)  and  (2-24b)  yields 


r(o)  a  / — —  exp[-i2Kn  (cose  )z]dz  (2-25) 

J  2nftcos V  0 

0  o  0 

? 

The  reflection  coefficient,  r(o),  times  cos  8Q  is  a  function  of 
cos  e_ 

— - — —  as  can  be  seen  from  Equation  (2-25).  There  are  several  places 
in  this  section  where  n(z)  is  assumed  to  be  of  the  form  given  by 
Equations  (1-1)  and  (1-2).  In  Appendix  A,  the  results  derived  in  this 
chapter  are  extended  to  n(z)  of  arbitrary  form. 

Now  divide  the  structure  shown  in  Figure  1-1  into  N  slabs  each 
parallel  to  the  x-y  plane  (see  Figure  2-1) 
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The  j  n  slab  has  thickness  d.,  and  refractive  index  n.(z)  given  by 

J  J 

j-1 


n j (z)  =  n(z  -  d£)  J  t 

£=  1 


n]  (z)  =  n(z) 


(2-26) 


Thus,  we  can  write 

0  2cos  6.  a 

2nQcos^  90  F(  - =  /n'(z)  exp[-i2Knn(cos  0jz]dz 


*/■' 


0  N  j-1 

■  ^ pj,xpC12  3  12-27a) 

j  =  2  £=1 


where 


and 


P,  =  I  nJ.(z)exp[-i2KnQ(cos  0Q)z]dz 
o 


♦j  *  Kn0(cos  0o)dj 


(2-27b) 


(2-27c) 


By  defining  n(z)  to  be  identically  zero  for  z>  D  and  z<  0,  we  can  write 


2cos  0, 


-f  * 


2nncos^  0„  F( - r — 2-)  =  /  n'(z)  exp[-i2Kn  (cos  0  )z]dz  (2-27d) 

U  O  A 


o'  o' 


2.3  TRANSFER  MATRIX  DEVELOPMENT 

Consider  a  pure  dielectric  material  which  has  a  refractive  index 
variation  along  the  z  axis  as  shown  in  Figure  2-2. 
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If,  as  indicated  in  the  Figure  2-2,  the  electric  fields  on  each  side 
of  the  dielectric  material  are  represented  as  the  sum  of  the  two  plane 
waves  propagating  in  opposite  directions,  then  EL  and  ER  can  be  related 
to  El  and  ER  through  a  two  by  two  transfer  matrix  S.  The  relationship 
is  given  explicitly  by 


(2-28) 


The  right-going  reflection  and  transmission  coefficients  r  and  t, 
respectively,  and  the  left-going  reflection  and  transmission  coef 
ficients,  r  and  t,  respectively,  are  defined  by 


r 


|r|exp[iY] 


ael 


(2-29a) 
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t  = 


r  = 


t  = 


A 

|t|exp[ir]  *  J- 

lR 


|r|exp[iy3  = 


|t|exp[ir]=  [± 


‘L®0 


ers° 


er  =  0. 


(2 -29b) 

(2-29c) 

(2-29d) 


If  follows  immediately  from  Equations  {2-28)  and  (2-29)  that 


r  = 

21 

(2-30a) 

t  * 

i 

Si 

(2-30b) 

r  = 

S12 

(2-30c) 

t  = 

C  S12S21 

22  “  S11 

(2-30d) 

sn 

-  1/t 

(2-31a) 

S12 

-  -  r/t 

(2-31 b) 

S21 

-  r/t 

(2-31c) 

CM 

CM 

CO 

=  t  -  rr/t 

(2-31d) 
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n 


Lt  tt  -  rrj 


Stoke’s  reversibility  theorem  requires  that 

t  r  +  r  t  =  0 


and 


t*t  +  rr*  =  1 


where  *  denotes  complex  conjugate. 

Equation  (2-33a)  implies 

r  =  -r*exp[i2r] 

Combining  Equations  (2-33b)  and  (2-34)  yields 

tt  -  rr  =  exp  [i2r] 

Finally,  with  Equations  (2-34)  and  (2-35),  Equation  (2-32) 


S  =  — 
a  t 


r*exp[i2r] 


exp[i2r ]J 


The  matrix  Sj  is  defined  by 


S  -  1 

SJ  *j 


rj 


rjexp[i2rjj 


exp[i2Tj]. 


(2-32) 


(2-33a) 

(2-33b) 


(2-34) 


(2-35) 

becomes 


(2-36) 


(2-37) 
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Sj  is  the  transfer  matrix  associated  with  the  j  dielectric  slab 
shown  in  Figure  2-1.  Thus,  the  transfer  matrix,  S,  of  the  entire 
structure  shown  in  Figure  2-2  is 


N 

n 

j=i 


sj  iri 

j=i  j 


where 


M  =  i 
hi  1 

C(M)  = 
L21 


rsj 

i 

shown 

that^] 

(N) 

r(N), 

11 

t12 

(N) 

r(N) 

21 

°22  J 

<  N 

(2-38) 


sum  of  all  op(m),  _<  N 


:{N)  =  $  exp[i2(|:  r  .)] 

1  J 


N 


:22}  =  Clf  exPti2<?  rj)3 


j=l 


and  where  ep(m)  is  an  expression  of  the  form 

m-1  *W*' 

VvVxpt1  s  2  0  -  <-l)j)r„] 

1  c  m  J-0  n=Jt. 

J 

A/rj  for  j  an  odd  integer 


bt  H 

j  Vj  for  j  an  even  integer 


(2-39) 


(2 -40a) 
(2 -40b) 
(2 -40c) 

(2-40d) 
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A 

r0  =  0 


m  is  a  positive  even  integer  and  1^,1^  •••» 

is  a  monoton i cal ly  increasing  sequence  of  positive  integers 


and 


op(m)  is  an  expression  of  the  form 


m-1 


£j+l_1 


>9  b  ...  b  exp[i  E  £  (1  - 

*1  j=0  n =£. 


(-i)j+1)rn3 


A  .r.  for  j  an  odd  integer 

b£  =  I 

j  V*  for  j 

J 


an  even  integer 

n„  4  o 

A 

ro  -  0 


(2-41 b) 


m  is  a  positive  odd  integer  and  2,^ ,  X^,  is  a  monotonically 

increasing  sequence  of  positive  integers. 

Using  Equations  (2-30a),  (2-30b),  (2-39)  and  (2-40),  we  can  now 


wri  te 


«(N) 


u2 1  _  sum  of  all  opCm),  ^m  -  N 
r(N)  "  1  +  sum  of  all  ep(m),x,m  ^  N 

(A  ,.) 


(2-42a) 


t  = 


'11 

i  _  _  ^j=i  y 

r(N)  1  +  sum  of  all  ep(m),  jl, 
C11  ™ 


r  _  sum  of  all  op(m) ,  Sn  1  N 


\M  V 


TK  <2'42b> 


(2-42c) 
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2erim 


As  an  example,  for  N=4  the  above  three  equations  become 


r^exp^rj)  +  r3exp(i2(r1  +  r2))  +  r4exp(i2(r1+  r2  +  r3)) 
+  r]rjr4exp(12  r3)  +  r2r*r4exp{i2(r1  +  r3))  +  r1r*r3exp(i2r2) 

+  r1r*r4exp(12(r2  +  r3)) 

1  +  r*r2exp(i2r1)  +  r3r4exp(i2r3)  +  r*r2r3r4exP(i2(ri  +  l*3)) 

+  r*r3exp{i2r2)  +  r*r3exp(i2(I^  +  r2)  +  r2r4exp(i2(r2  + 

+  r*r4exp(i2(r1  +r2  +r3)) 


t2t3t4exp ( i  ( +  ^2  +  ^3  +  ^4^ 

1  +  r*r2exp(i2r1)  +  r*r4exp(i2r3)  +  r*r2rjr4exp(12(rj  +  r3)) 
+  r*r3exp(i2r2)  +  r*r3exp(i2(T1  +  r2))  +  rjr4exp(i2(r2  +  r3)) 
+  r*r4  exp(i2(r1  +  r2  +  r3)) 


r-j  +  r2exp(i2r1)  +  r3exp(i2(r1  +  r2) )  +  r4exp(i2(r1  +  r2  +  r3)) 

+  r1r*r4exp(i2r3)  +  r2r*r4exp(i2(r1  +r3))  +  r^^exp^y 

r  +  r1r2r4exp^i2^r2  +  r3^ 

*  tlt2t3t4 

Formulas  analogous  to  Equations  (2-42a)  and  (2-42b)  were  first  derived 
for  multilayer  dielectric  stacks  b,  Crook^18^  Equations  (2-42)  have 
been  exploited,  by  Pegis^8^,  Delano^"8^,  and  Knittl^  to  aid  in  the 
design  of  mulitlayer  dielectric  stacks  and  inhomogenous  layers. 
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2.4  THE  BASIC  DESIGN  EQUATION 

When  n(z)  is  given  by  Equations  (1-1)  and  (1-2)  and  when 

p 

\r.(z0)\  «  1  for  0  <  zQ<  d^,  then  from  Equation  (2-25) 


and 


-j  ■/ 


(z) 


ocos  60 


2 —  exp[-i2Kn0(cos  eo)z]dz  (2-43) 


f  ^ 

Jd  2n  co; 


-  z) 


/dj  2nocos  eo 


exp[-i2Kn0(cos  eQ) (d^  -  z)]dz 


=  -exp[-i2KnQ(cos  e^] 


2n  cos  e 
o  o 


2 —  exp  [i2KnQ(cos  60)zjdz 


“  -  r*exp[-i2KnQ(cos  e^.]  (2-44) 

Combining  Equations  (2-43),  (2-44)  and  (2-34)  yields 


Tj  =  -iKnQ(cos  e0)dj  (2-45) 

when  |rj(zQ)|2«l,  0<zo<d... 

2 

We  note  that  the  condition  Jr^. (zQ)  |  «  1,  0  <_  zQ  <_  d^  can  always  be 

satisfied  by  choosing  d.  sufficiently  small. 

J 

According  to  (2-41 b), 

(sum  of  all  op ( 1 ) ,  ^  <:  N)=  r-j  +  r2exp[i2r-|]  +  r3exp[i2(r,  +  r2)] 

+  ...  +  rNexp[i2(r1  +  r2  +  ...  r^)] 
(2-46) 
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Now  from  Equations  (2-43),  (2-45),  (2-46),  and  (2-27),  it  follows  that 

00 

2n0(cos2e0)(sum  of  all  op(l),£1£  N )-J  n‘ (z)exp[-2KnQ(coseo)z]dz 


(2-47) 


when  I  rj(z0)r  «  0  <  z0  1  dj 


According  to  Equation  (2-41),  op(m)  is  a  product  involving  an 
odd  number  (m)  of  r ^ ’ s .  Similarly,  ep(m)  is  a  product  involving  an 
even  number  (m)  of  r.'s.  When  the  d.'s  are  small,  then  the  r.'s  will 
be  small.  Thus,  we  expect  |op(m)|  and  |ep(m)|  to  decrease  rapidly  as 
increases.  If  all  op(m)  and  ep(m),  m  >  1  are  neglected  then  Equation 
(2-42a)  becomes 


r  =  sum  of  all  op(l),  Jtj  £  N  (2-48a) 

If  all  op(m)  and  ep(m),  m  >  2  are  neglected  then  Equations  (2-42a)  and 
(2-42b)  become 


„  sum  of  all  op(l),  &1  -  N 


(2 -48b) 


r  __  sum  of  all  op(1 ) ,  &1  -  N 


(2 -48c) 


It  is  expected  that  Equations  (2-48b)  and  (2-48c)  are  more  accurate 
than  Equation  (2-48a)  since  fewer  terms  were  neglected  in  their 
derivation. 

Recall  that  n(z)  is  linked  to  op(l)  through  Equation  (2-47).  Thus, 
Equations  (2-48a)  and  (2-48c)  can  be  rewritten  as 
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oc 

./• 


2n0(c°s  0Q)r  ~ /  n’(z)exp[-i2Kn0(cos  0Q)z]dz  (2-49a) 


and 


2nQ(cos  e 


'(z)exp[-i2Kno(cos  0Q)z]dz  (2-49b) 


respectively.  Equation  (2-48b)  appears  to  be  of  little  use,  since  a 
simple  expression  for  the  sum  of  all  op(2),  <  N  in  terms  of  n(z) 

is  not  known.  The  right  hand  sides  of  Equations  (2-49a)  and  (2-49b) 
are  the  Fourier  transform  of  n'(z).  By  the  property  of  Fourier  transforms 


j  n'(z)exp[-i2Kn0(cos  0Q)z]  - 

.00 


i2nQ(cos  0Q) 
X 


z)exp[-i2KnQ(cos  0Q)z]dz 


(2-50) 


Using  Equation  (2-50),  noting  that  n(z)  is  a  real  function,  and  taking 
the  inverse  Fourier  transform  of  both  sides  of  Equations  (2-49a)  ar.d 
(2-49b),  yields 

oo 

/i  i  2n  cos  0 

c-0-s-  q  r  exp[-2KnQ(cos  e0)z]d( - - — 5.)  (2-51a) 

o  0 


and 

n(z)  “  2 Re 


°°  N 

/-lx  r  (IT  tX  exp[-iKn  (cos  en)z]d( 
cos"e7t\j=l  J/  0  0 


2nQcos 


-i) 

(2-51 b) 


respectively.  In  Equations  (2-51  a)  and  (2-51b),  Re  denotes  the  real 

part.  Equation  (2-51a)  is  only  valid  in  the  regime  of  low  reflectivity, 
o 

i.e.  small  |r|  ,  whereas  Equation  (2-51b)should  also  be  valid  in  regimes  where 
the  reflectivity  is  not  low.  Equation  (2-51 b)  is  the  central  result  of 
this  chapter  and  will  form  the  basis  of  our  iterative  Fourier  Transform 
reflection  grating  design  technique. 
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The  Fourier  Transform  relation  between  r  and  n(z)  for  small  r, 
given  by  Equation  (2-51  a )  is  well  known^20^. 


A  2 

In  our  design  technique,  R  =  I r I  ,  is  assumed  to  be  given  as 
cos  @o 

a  function  of - r — ’  and  the  objective  is  to  find  the  refractive 

A  COS  0 

index,  n(z),  which  results  in  this  R  vs.  — ^-characteristic. 
Since  n  =  n  (see  Figure  (1-l))energy  conservation  dictates  that 


jr|2  +  |tj2  =  1 


(2-52) 


Thus,  Equation  (2- 51b)  can  be  rewritten  as 


,  R  1/2  n 

„(2)  .  2Re  /  -  nP-e 


/A"  \  2n  cos  6 

a(jUl  |tj  |yexp[i2Knncos  ejz]d(— 5- 


o - o' 

(2-53) 


where 


a  =  exp[i(y  -  r  +  E  r  )]  (See  Equations  (2-29))  (2-54a) 

Jl=l 


=  exp[i(Y  -  r  -  KnQ(cos  0Q)D]  (See  Equation  (2-45)  and  Figure  (2-1)) 

(2-54b) 


Y  and  r  are  the  phases  of  the  reflection  (r)  and  transmission  (t) 

coefficients,  respectively.  Also  note  that  a  is  a  function  of 
cos  e  cos  0 

- -  .  In  general,  R  vs.  — * — -  will  be  specified  by  the  design, 

A  cos  9  A 

and  a  vs.  - will  not  be  specified. 

X 

Our  design  procedure  goes  as  follows: 


1. 

2. 

p. 

Iteration 
Process  4 

5. 

-6. 


A  refractive  index,  n(z),  is  chosen. 

a  is  assumed  to  be  identically  1. 

Using  Abele's  multilayer  theory,  the  t/s  are  computed 
from  n(z).  J 

n(z)  is  calculated  from  Equation  (2-53). 

The  corresponding  R  vs.  cos  eQ/A  characteristic  Is  computed. 
Step  3  is  repeated  until  nui  changes  little  from  one 
iteration  to  the  next. 
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7.  The  resulting  n(z)  is  the  desired  design. 

2.5  DESIGN  EXAMPLES 

The  procedure  outlined  above  was  used  to  design  two  reflection 
grating  filters.  The  first  filter  was  specified  as  having  the 
reflectance  vs.  wavelength  characteristics  shown  in  Figure  [2-3a], 
while  the  second  filter  was  specified  as  having  the  characteristic 
shown  in  Figure  [2- 3b] .  The  evolution  of  the  designs,  as  a  function 
of  the  number  of  iterations,  is  shown  in  Figures  [2-4  and  2-6].  Given 
in  Figure  [2-5]  is  the  refractive  index  profile  corresponding  to  the 
last  iteration  shown  in  Figure  [2-4].  Similarly  Figure  [2-7]  gives  the 
refractive  index  profile  corresponding  to  the  last  iteration  shown  in 
Figure  [2-6].  The  refractive  index  profile  is  displayed  by  plotting 
both  An(z)  (solid  curve)  and  f(z)  (dashed  curve)  as  a  function  of  z. 

For  both  designs,  nQ  =  1.5,  grating  thickness  =  180  ^m,  and  en  =  0°. 

The  initial  choice  of  n(z)  for  both  filters  one  and  two  was 

n(z)  =  1.5  +  0.015  sin  2iT(5.88um-^ )z 

An  examination  of  the  points  plotted  in  Figures  [2-4]  and  [2-5]  would 
indicate  that  both  f(z)  and  An(z)  are  essentially  constant  over  intervals 
of  length  l/f(z).  Thus,  n(z)  Is  simply  a  sinusoid  whose  frequency  and 
amplitude  are  slowly  changing  as  a  function  of  z.  The  frequency  "spikes" 
shown  in  Figure  [2-7]  occur  at  points  where  An(z)  =  0.  Thus,  these  spikes 
correspond  to  spacers  (i.e.,  regions  where  An(z)  =  0)  of  length  l/(fz). 


Although  a  detailed  description  of  the  computer  implementation  of 
Equation  (2-53)  is  not  provided  in  this  report, some  important  points 


FIGURE  2-5.  REFRACTIVE  INDEX  PROFILE  OF  A  GRATING  HAVING  A  TRIANGULAR 
CHARACTERISTIC.  Solid  Curve  is  An(z)  and  Dashed  Curve  is  f(z). 


20.00 


FIGURE  2-7.  REFRACTIVE  INDEX  PROFILE  OF  A  GRATING  HAVING  A  SADDLE 
CHARACTERISTIC.  Solid  Curve  is  An(t)  and  Dashed  Curve  is  f(z). 
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associated  with  this  implementation  will  be  discussed.  First,  it 
is  noted  that  the  n(z)  obtained  from  Equation  (2-53)  will  be  of 
infinite  extent  in  z.  After  each  iteration,  one  must  truncate  the 
extent  of  n (z) ,  so  that  An(z)  is  nonzero  only  over  the  interval  of 
finite  length  corresponding  to  the  filter  thickness.  Generally 
An(z)  is  quite  small  for  large  values  of  z,  and  thus,  little  accuracy 
is  lost  by  the  truncation.  Second,  in  order  to  compute-^  |t.|» 

the  thickness,  d.  of  each  slab  must  be  chosen.  The  approximations 

J 

used  in  the  derivation  Equation  (2-53)  are  sensitive  to  these  d's. 

J 

In  our  designs,  for  each  value  of  cos  eQ/x  at  which 


(-hj 


1/2  N 


j«l 


p 

was  evaluated,  di  was  chose*  so  that  |t,|  =  0.075.  Thus,  the  d.'s  are 

J  J  J 

not  fixed  during  the  design  of  a  given  filter,  but  change  for  each  value 
of  cos  0qA.  The  optimum  choice  for  the  dj's  has  not  been  determined. 
Third,  the  t.'s  associated  with  each  iteration  and  the  resulting  R  of 

J 

the  final  design,  are  evaluated  using  an  efficient  computer  imple¬ 
mentation  of  Abeles'  multilayer  matrix  technique.  Details  of  this 
implementation  can  be  found  elsewhere^7^.  At  each  iteration  an  n(z) 
is  known,  and  so,  Abells'  multilayer  matrix  technique  can  be  used  to 
determine  y  and  r,  and  hence  a  .  Thus,  it  is  noted  that  a  better 
design  would  most  likely  be  obtained  by  computing  a  at  each  itera¬ 
tion  rather  than  assuming  that  a  is  identically  one. 
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APPENDIX/SECTION  2 


In  this  appendix,  the  reflection  grating  design  technique  is 
extended  so  as  to  be  applicable  to  multilayer  dielectric  stack 
design.  The  multilayer  dielectric  stack  design  problem  consists 
of  two  phases.  First,  the  refractive  index  profile,  n(z),  required 
to  realize  a  specified  reflectance  vs.  wavelength  characteristic  at 
a  fixed  angle  of  incidence,  is  determined.  Second,  a  discontinuous 
profile,  which  could  be  physically  realized  by  a  multilayer  dielec¬ 
tric  stack,  is  fitted  to  the  continuous  refractive  index  profile, 
n(z).  The  discussion  here  is  restricted  to  the  first  phase  of  the 
above  mentioned  design  problem.  The  technique  for  determining  n(z) 
given  the  reflectance  vs.  cos  6Q/\  has  already  been  developed  in 
section  two.  This  technique,  however,  required  n(z)  to  be  of  the  form 
given  by  Equations  (1-1)  and  (1-2).  We  will  remove  this  restriction 
below  and  allow  for  arbitrary  n(z).  The  extension  to  arbitrary  n(z) 
involves  only  a  small  modification  of  our  reflection  grating  design 
technique  and  so  our  discussion  will  be  brief. 


Repeated  below  is  Equation  (2-18)  which  it  will  be  recalled  is 
valid  for  arbitrary  n(z)  provided  that  |r(z)  j  «  1 ,  0  <_  z  <_  D. 


where 


r(o) 


y(z) 


D  P 

=/fe^yexpC'i2K 

O  0 

=  n(z)cos  e(z) 


(2-18) 

(2-21 b) 


eU) 


n(z)cos  e(z) 


( 2 -2 1  c ) 
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Next,  n  is  defined  by 


Z 

n  =  2fy(i)dz  (A-l ) 

o 

It  is  noted  thatn(z0)  is  twice  the  optical  path  length  that  a  plane 
wave  travels  into  the  dielectric  in  order  to  reach  the  plane  z  =  zQ. 
Combining  Equations  (2-18),  ( 2 -21 b)  and  (2-21 c )  yields 


r(o): 


/ 


-n  ^ ■  -  ■  exp[-iKn]dn 


2rT(Ti)cos  "Stn) 


(A-2) 


where 


0 

f 


L  =  2  /  Y(  t>dc 
0 


(A-3) 


n(n)  =  n(z) 


and 


¥(n)  =  6(z) 


n  and  z  are  related  via  Equation  (A-l) 

(A-4) 


Thus,  fT(n)  and  ¥(n)  are  simply  the  refractive  index  and  propagation 
angle,  respectively,  expressed  as  a  function  of  twice  the  optical 
path. 

C(n)  and  F  are  defined  by 


C(n ) 


n(n) _ 

-  2  - 
2n(n )  cos  e(n) 


_ n(n)n,(n) _ 

Z(n2(n)  -n2(0)  sin2  0(0)) 


L 

C(n)exp[-iKn]dn 


(A-5) 


(A-6) 
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It  follows  immediately  from  Equations  (A-5)  and  ( A-6 )  that 


where 


N 

M 

JVexp[i 

Eu 

(A-7) 

j=2 

m=l 

N 

Zd.  =  L 
j  =  l  3 

(A-8) 

A  / 

P.  =  /  C(  n)exp[-iKri]dn 

(A-7) 

*i  '  Kd3 

(A-10) 

Using  Equation  (A-2)  it  is  easy  to  show  (see  for  instance  Equations 
(2-43)  and  (2-44))  that 


i"j  =  -  r^.expC-iKdj]  (A-ll) 

Combining  Equations  (A-ll)  and  (2-34)  now  yields 


rj*~iKdj/2  (A-12) 

From  Equations  (2-46),  (A-12)  and  (A-7)  -  (A-9) 

F  =  sum  of  all  op(l),  t-|  ( A-l 3) 

From  Equations  (A-l 3) ,  (A-6)  and  (2-48c)  it  follows  that 


,  N  « 

(‘  P  **):/ 


C(n)exp[-iKnJdn 


(A-l 4) 
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Where  C(n)  is  defined  to  be  identically  zero  forn<  0  orn  >  L. 
Delano^"®-!  derives  an  expression  similar  to  Equation  (A-14),  but 
assumes  as  a  first  order  approximation  that 


N 

TT 

j=i 


t,  S  1 


From  Equation  (A-14) 
C(n)  =  2Re 


00  ^ 

AfevexpMMd(‘) 


(A-15) 


where  r  and  t  are  specified  as  a  function  of  1/A  and  the  t.'s  are 

J 

computed  as  a  function  of  1/*.  Note  that  Equation  (A-15)  is  analogous 
to  Equation  (2-51 b) .  Once  C(n)  is  computed  from  Equation  (A-15), 
the  differential  equation  defined  by  (A-5)  can  be  solved,  numerically, 
to  obtain  n(n). 


This  differential  equation  can  be  solved  analytically  for  the  case 
of  normal  incidence,  i.e.  9(o)  *=  0.  Finally,  rf  is  obtained  from  n  by 
using  the  following  relationship 


(/U_). 

yo  nU)cos  ¥(?)/ 


n(n) 


(A- 16) 


or  equivalently  (see  Equation  (2-6a)) 

dc 


'(/i 


[^(0  -  nZ(o)sin20(o)]V/2 


(n)  (A-l  7) 


We  note  that  Equation  (A-15)  can  be  rewritten  as 


oo 

/1/2 

(f-;  fi)  a  *TT|tj  |exp(iKn)d(-|-) 


(A-l 8) 
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where 


N 

a  =  exp[i (Y  -  r  +  r-)3  (A-19) 

j=l  J 


Equation  (A-18)  is  analogous  to  Equation  (2-53). 

Delano  distinguishes  between  two  types  of  synthesis  problems, 
ones  where  R  is  given  (first  type),  and  ones  where  r/t  is  given 
(second  type).  He  handles  these  two  types  of  problems  differently.  We 
handle  synthesis  problems  of  the  first  type  via  Equation  (A-15)  and  those 
of  the  second  type  via  Equation  (A-18).  Thus,  it  is  seen  that  our  tech¬ 
nique  draws  little  distinction  between  the  two  types  of  synthesis  prob¬ 
lems.  Sossi  [9]  derives  an  equation  similar  to  Equation  (A-18)  for  the 
case  of  normal  incidence  (i.e.,  0Q  »  0).  In  his  equation  the  quantity 


1/2 


We  have  been  unable  to  follow 


tj  is  replaced  by  (1  -  £  ) 
his  derivation.  Implicit  in  our  derivation,  in  this  appendix,  has  been 
the  assumption  that  the  refractive  index,  n(z),  does  not  have  disconti¬ 
nuities  at  the  boundaries  z  =  0  and  z  =  D  (i.e.,  nc  =  n(o)  and  n$  = 
n(D)).  We  believe  that  our  technique  can  be  modified  to  remove  this  re¬ 
striction.  It  is  also  noted  that  our  technique  can  be  easily  extended 
to  handle  the  case  of  p,  rather  than  s,  polarization. 
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Section  3 


CORRUGATED  WAVEGUIDE  FILTERS 


A  formal  mathematical  analogy  between  reflection 
gratings  and  corrugated  waveguide  filters  (CWF)  is 
demonstrated  in  this  chapter.  The  possibility  of  de¬ 
signing  CWF,  using  the  iterative  Fourier  transform 
technique  developed  in  Section  2,  is  explored. 
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3-1.  INTRODUCTION 

In  this  section  we  demonstrate  a  formal  mathematical  analogy  between 
reflection  gratings  and  corrugated  waveguide  filters  (CWF).  The  devel¬ 
opment  given  by  Yariv  [1]  is  closely  followed.  We  emphasize  that  the 
mathematics  developed  for  CWF  are  not  rigorous.  Consequently,  predicted 
results  should  be  experimentally  verified  and  attempts  should  be  made  to 
develop  a  rigorous  mathematical  approach. 

3-2.  FORMAL  MATHEMATICAL  ANALOGY  BETWEEN  REFLECTION  GRATINGS  AND  CWF 

Consider  the  CWF  shown  in  Figure  [3-1].  The  filter  consists  of  a 

thin  waveguiding  layer  (0  <  x  <  a)  of  refractive  index  ng  sandwiched 

between  two  dielectrics  having  refractive  indices  n  and  n  (n  >  n 

c  s  g  s 

>  nc).  The  upper  surface  of  the  waveguiding  layer  has  a  corrugation  of 
small,  maximum  depth  (i.e.,  hQ  «  a).  The  refractive  index,  n(x,  z), 
can  be  written  as 

n2(x,  z)  =  n2  +  w(x ,  z)(n2  -  n2)  sin  27rf(z)z  (3-1) 

where  nb  equals  nc  in  region  C,  n^  in  region  G  and  n$  in  region  S, 
and  wherew(x,  z)  is  a  binary  function  which  assumes  either  the  value 
zero  or  one.  Note  that  w(x,  z)  =  0  for 

x  <  0  or  x  >  a 

If  we  restrict  our  discussion  to  transverse  electric  fields,  then  Max¬ 
well's  wave  equation  in  all  three  regions  can  be  written  as 

2 

V2Ey(x,  z)  +  (^)  n2(x,  z)Ey(x,  z)  *  0  (3-2a) 

or  equivalently 
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2  2 

v2Ey(x,  z)  +  0)  n2Ey(x,  z)  +  0)  w(x,  z)(n2  -  n2)  sin  (27rf(z)z)  =  0 

(3-2b) 

Since  the  corrugation  depth  is  small,  we  will  assume  that  the  total  field 
in  the  "perturbed"  waveguide  can  be  approximately  written  as  the  super¬ 
position  of  the  confined  modes  of  the  waveguide  without  the  corrugation, 
i.e., 

Ey(x,  z)  =  £Ym)(z)E^m)(x)  exp  [-ip[m)z] 
m 


+  B(m)(z)E<m)(x)  exp  [iP[m)z] 


(3-3) 


In  Eq.  (3-3),  Eym^(x)  exp  C+ip confined  mode  of  the  un¬ 
corrugated  waveguide.  Thus 


(m). 


th 


1  E^'lx)  exp  +(2jL)2„2Eh)(x)  exp  [♦1P<">2] 


=  0 


(3-4) 


and  is  the  propagation  constant  of  the  mth  mode  and  can  be  written  as 


(m)  _  2tt_ 


,  (m)  - 


'Z 


(3-5) 


where  e^(x)  is  the  propagation  angle  of  the  mth  mode  within  the  guiding 

layer.  Note  that  e£m),  P^m)  and  E^m)  are  all  functions  of  x.  Given  the 

geometry  of  the  uncorrugated  waveguide,  0^,  p^  and  are  easily 
determined  [1].  For  the  sake  of  simplicity  the  waveguide  will  be  assumed 
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to  be  single  mode.  Equation  (3-3)  now  reduces  to 

Ey(x,  z)  =  A^(z)E^(x)  exp  C-ip^1^]  +  ^(z)Ey^  ^(x)  exp  [ip^z] 

(3-6) 

The  reflection  and  transmission  coefficients  r  and  t  of  the  CWF 
are  defined  as  B^^(0)/A^\o)  and  A^1  ^(O/A^1  ^(o)  respectively.  A 
simple  physical  argument  indicates  that  B^(L)  =  0.  Substituting  Eq. 
(3-6)  into  Eq.  (3-2b),  and  then  using  Eq.  (3-4),  yields 

A(1,"(z)e|1)(x)  exp  [-ip^z]  -  21P<1)A(1),(z)E<1)(x)  exp  [-ip^z] 

+  B^*  (z)E^(x)  exp  [ip^z]  +  2ip^  B  ^  (z)E^(x)  exp  [ipj^z] 

=  -(~)  w(x,  z) (n^  -  n^)  sin  (z-nf(z)z)  [V1  ^(z)E^(x)  exp  (-ip^z) 

+  B(1)(z)E^1}(x)  exp  (ip^zjj  (3-7) 

Now  it  will  be  assumed  that  Ey^(x)  has  been  normalized  so  that 


J  [^(x)]2  dx  =  1 


(3-8) 


Furthermore,  since  the  maximum  corrugation  depth,  hQ,  is  small  compared 
to  the  thickness,  a,  of  the  waveguiding  layer 


5 


1 


RADAR  AND  OPTICS  DIVISION 

E^1}(X)  =  E^r1)(a)  for  w(x,  z)  f  0  (3-9) 

Multiplying  both  sides  of  Eq.  (3-7)  by  E^(x)  and  then  integrating  from 
to  «  yields 

A^  (z)  exp  [-ip^z]  -  2ip^(z)A^  (z)  exp  [ip^z] 

+  (z)  exp  [ip^z]  +  2ip^V^  (z)  exp  [ip^z] 

2  2  °° 

-  -(t)  (nc  '  "g>A(1)<z)  sin  (2uf(z)z)  exp  (-ip*1  }z) TeJ1  }(a)]  J  w(x,  z)  dx 

*  **  -OD  * 

00 

-  (t)Z{\  -  ng>B(1)(Z>  Sin  *2"f*Z)Z)  exP  ^M2I  W(X>  Z)  dX 

*00 

(3-10) 


where  we  have  used  Eqs.  (3-8)  and  (3-9).  Recall  that  E^(a)  is  a  func¬ 
tion  of  x.  The  functions  e(x)  and  h(z)  are  defined  by 


(3-11) 
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CO 

=  j  w(x,  z)  dz 


(3-12) 


Observe  that  h(z)  is  the  corrugation  height.  Using  the  above  two 
definitions.  Eg.  (3-10)  becomes 


A^  (z)  exp  [-ip^z]  -  2ip^A^  (z)  exp  (-ip^z] 


+  (z)  exp  [ip^z]  +  2ip^1V1^  (z)  exp  [ip^zj 


2  2 


-  exp  [1(2.f( 


x)  -  P^’lx] 


_2  2 


+i^)  ng~ e(x)h(z)A^(z)  exp  [-i(2irf(z)  +  p^)z] 


2  n2  2 

’(jf)  ng  - - 1  e(x}h(z)B^(z)  exp  [i(2irf(z)  +  pj^z] 


/  \2  r?  -  n^ 

i(f-)  n  c  ---  9  e(x)h(z)B(1)(z)  exp  [-1(2irf(z)  -  P<1})z]  *  0  (3-13) 

W  9  2n_  z 


Now  Eq.  (3-13)  is  identical  to  Eq.  (1-8)  if  the  following  quantities 
are  equated. 
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> 


Thum 


A(z)  ~  A(1)(z) 

(3-14a) 

B(z)  B(1 *(z) 

(3-14b) 

n0~  ng 

(3-1 4c ) 

e0~  eg1)(x) 

(3-14d) 

pz  +-+  pz^  (x) 

(3-1 4e ) 

n2  -  n2 

An ( z )  g  e(x)h(z) 

( 3-1 4f ) 

2"g 


Recall  that  for  any  given  X ,  e*"(x),  e(x)  and  o*’*(x)  are  known.  Thus, 
if  h(z)  is  also  known,  Abelfes'  multilayer  dielectric  theory,  as  pre¬ 
sented  in  Chapter  can  be  used  to  find  r.  Let  us  divide  the  CWF  into 

N  sections  along  the  z  axis  with  the  reflection  and  transmission  coef- 
h 

ficients  of  the  j  section  being  denoted  by  r.  and  t.,  respectively. 

J  J 

Then,  as  in  Section  3, 

N 

tll^  =op(1)  (3-i5) 


and 


op(l ) 


•/ 


2  2 
n_  -  n 


2 — - e(x)h'(z)  exp  [i2Kn  (cos  ej^(x))z]  dz 

c  qV  1  >f\  \  9  9 


2n"  cos"  ey  ;(x) 
9  9 


(3-16) 


I 


i 

i 

i 
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Combining  Eqs.  (3-15)  and  (3-16)  yields 


h'  (z) 


^  cos2  .»>|» 

i-i  /  ("c  -  $«<*> 


2n„  cos  e^(> )' 


=  /  tflltj)  6XP  t^Kng(cos  ej1>(x))x]d(-a— J 


(3-17) 


Using  the  properties  of  Fourier  transforms  Eq.  (3-17)  becomes 
N 


h(z)  =  2Re 


r  /  "  \  2in  cos  e^(x) 

I  I  t .  I  — - exp 

{  'j-i  7  (nc  -  4) e(i) 


/  n  \  \  /<-n„  cos  el^(A) 

[i2Kng(cos  eg  (x)^z]d' — ^ ^ 


Thus,  if  r  and  t  are  known  as  a  function  of  X,  then  the  iterative  tech¬ 
nique  outlined  in  Section  3  can  be  used  to  determine  h(z). 
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